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Abstract 

Jet cross sections were measured in charged current deep inelastic e^p scatter- 
ing at high boson virtualities with the ZEUS detector at HERA II using an 
integrated luminosity of 0.36 fb~^. Differential cross sections are presented for 
inclusive-jet production as functions of Q^, Bjorken x and the jet transverse en- 
ergy and pseudorapidity. The dijet invariant mass cross section is also presented. 
Observation of three- and four-jet events in charged-current e^p processes is re- 
ported for the first time. The predictions of next-to-leading-order (NLO) QCD 
calculations are compared to the measurements. The measured inclusive-jet cross 
sections are well described in shape and normalization by the NLO predictions. 
The data have the potential to constrain the u and d valence quark distributions 
in the proton if included as input to global fits. 
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1 Introduction 



Jet production in charged- current (CC) deep inelastic e^p scattering (DIS) provides a 
testing ground for QCD and for the electroweak sector of the Standard Model (SM). Up 
to leading order (LO) in the strong coupling constant, ctg, jet production in CC DIS 
proceeds via the quark-parton model {Wq — > g, Fig. [1^), ly-gluon fusion {Wg qq, 
Fig. db) and the QCD-Compton {Wq — > qg, Fig. [lb) processes. Thus, differential cross 
sections for jet production are sensitive to both the value of and the mass of the 
propagator, Mw, which are fundamental parameters of the theory. Cross sections in CC 
DIS are also sensitive to the valence flavor content of the proton, since the W~^~^^ couples 
to the up-like (down-like) quarks in the proton. 

The large center-of-mass energy available at the HERA e^p collider {^/s = 318 GeV) 
extends the kinematic region for studying CC DIS with respect to fixed-target neutrino 
scattering experiments [1] by about two orders of magnitude in the virtuality of the 
exchanged boson, Q^, and to lower values of the fraction of the proton momentum carried 
by the struck parton, x. Measurements of the CC DIS cross section at HERA [2-4] 
demonstrated, at high Q^, the presence of a space-like propagator with a finite mass, 
consistent with that of the W boson. During 2002-2007, HERA provided longitudinally- 
polarized electron or positron beams. Measurements of the fully-inclusive CC DIS cross 
section for positive and negative values of the longitudinal polarization of the beams were 
found to be in good agreement with the predictions of the SM [5,6]. 

At HERA, multijet structures were observed in CC DIS [3,7,8] at large and jet shapes 
and subjet multiplicities were measured [8,9] and compared with neutral current (NC) 
DIS processes [9,10]. The subjet multiplicities were used to extract a value of as(M^) [8]. 

In this paper, measurements are presented of inclusive-jet and dijet cross sections in CC 
e^p DIS in the laboratory frame. Measurements of three-jet differential cross sections in 
CC DIS were measured for the first time in e^p collisions. A small sample of four-jet 
events was also observed in the data. The measurements are presented as functions of 
Q^, X, the jet transverse energy, and pseudorapiditjQ, rj^'^^, and as a function of the 
invariant mass of the jet system in dijet and three-jet events. Predictions from next- 
to-leading-order (NLO) QCD calculations are compared to the measurements. Results 
for negative and positive longitudinally-polarized electron and positron beams are also 
presented. 



The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the center of HERA. The coordinate origin is at the nominal interaction point. 
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2 Experimental set-up 



A detailed description of the ZEUS detector can be found elsewhere [11,12]. A brief 
outline of the components most relevant for this analysis is given below. 

Charged particles were tracked in the central tracking detector (CTD) [13], which operated 
in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The CTD 
consisted of 72 cylindrical drift-chamber layers, organized in nine superlayers covering the 
polar-angle region 15° < 6 < 164°. In 2001, a silicon microvertex detector (MVD) [14] was 
installed between the beampipe and the inner radius of the CTD. The MVD was organized 
into a barrel with three cylindrical layers and a forward section with four planar layers 
perpendicular to the HERA beam direction. The barrel contained 600 single-sided silicon 
strip sensors each having 512 strips of width 120 /im; the forward section contained 112 
sensors each of which had 480 strips of width 120 /im. Charged-particle tracks were 
reconstructed online by using the ZEUS global tracking trigger [15], which combined 
information from the CTD and MVD. The online tracks were used to reconstruct the 
interaction vertex and reject non-ep background. Offline, the tracks used in this analysis 
were reconstructed using information from the CTD and were used, in addition to the 
vertex reconstruction, to cross-check the energy scale of the calorimeter. 

The high-resolution uranium-scintillator calorimeter (CAL) [16] covered 99.7% of the 
total solid angle and consisted of three parts: the forward (FCAL), the barrel (BCAL) 
and the rear (RCAL) calorimeters. Each part was subdivided transversely into towers 
and longitudinally into one electromagnetic section (EMC) and either one (in RCAL) 
or two (in BCAL and FCAL) hadronic sections (HAC). The smallest subdivision of the 
calorimeter is called a cell. Under test-beam conditions, the CAL single-particle relative 
energy resolutions were a{E)/E = 0.18/V^ for leptons and a{E)/E = 0.35/V^ for 
hadrons, with E in GeV. 

The luminosity was measured using the Bethe-Heitler reaction ep —>■ e^p by the luminosity 
detector which consisted of two independent systems. In the first system, the photons 
were detected by a lead-scintillator calorimeter placed in the HERA tunnel 107 m from 
the interaction point in the lepton-beam direction. The system used in previous ZEUS 
publications [17] was modified by the addition of active filters in order to suppress the 
increased synchrotron radiation background of the upgraded HERA collider. The second 
system was a magnetic spectrometer arrangement [18], which measured electron-positron 
pairs from converted photons. The fractional uncertainty on the measured luminosity was 
3.5%. 

The lepton beam in HERA became transversely polarized naturally through the Sokolov- 
Ternov effect [19], with a build-up time of approximately 40 minutes. Spin rotators 
on either side of the ZEUS detector changed the transverse polarization of the beam 
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into longitudinal polarization. The lepton-beam polarization was measured using two 
independent polarimeters, the transverse polarimeter (TPOL) [20,21] and the longitudinal 
polarimeter (LPOL) [21,22]. Both devices exploited the spin-dependent cross section 
for Compton scattering of circularly polarized photons off leptons to measure the beam 
polarization. The fractional uncertainty on the measured polarization was 4.2% and 3.6% 
from TPOL and LPOL, respectively. 



3 Data selection and jet search 

The data were collected from 2004 to 2007, when HERA operated with protons of energy 
Ep = 920 GeV and electrons (positrons) of energy = 27.5 GeV, and correspond 
to an integrated luminosity of 180.0 ± 6.3 (178.5 ± 6.2) pb"^. Samples of negatively- 
(positively-) polarized electron beams with an integrated luminosity of 106.4 (73.6) pb^^ 
and luminosity-weighted average polarization of P^^^^^ = — 0.27±0.01 {P^-^ = +0.29±0.01) 
were analyzed. For positrons, the samples analyzed were of 76.5 and 102.1 pb~^ with a 
luminosity-weighted average polarization of P^+^ = — 0.371o;q2 and P^°^ = +0.32 ± 0.01, 
respectively. 

The main signature of a CC DIS event at HERA is the presence of large p^^^^ and large 
E^^, where p^^^^ is the missing transverse momentum arising from the energetic final-state 
neutrino which escapes detection and E^^ is the total transverse energy arising from the 
hard interaction. The variable p™'^^ was reconstructed using the vectorial sum p™'**** = 
+ Cl2iPY,iy and E^^ was reconstructed as the scalar sum E^^^ = Et^i- In 
both cases, the sum runs over all GAL cells. The online selection of the signal was based on 
the ZEUS three-level trigger [12]. Two different trigger selections were used. One trigger 
selection relied on large p™*^*^ and large E^^. The alternative trigger selection additionally 
required the presence of at least one jet with transverse energy above 8 GeV. Events 
from GG DIS interactions were selected offline using criteria similar to those of an earlier 
publication [8]. The kinematic variables Q^, the inelasticity, and x were estimated 
using the method of Jacquet-Blondel [23], which uses the information from the hadronic 
energy flow of the event, and corrected for detector effects as described elsewhere [9,24]. 
These estimators were reconstructed as: 

ym = t;-^ , Qm = and xjb - 



2E, 1-yjB ' syjB 

where the sum runs over all GAL cells. The main selection criteria were: 

• pf'' > 11 GeV; 

• p'^^^^/E^^ > 0.5, to reject photoproduction and beam-gas background; 
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• the vertex position along the beam axis in the range —35 < Z^t^ < 33 cm, consistent 
with an e^p interaction; 

• at least one track associated with the vertex, which had a polar angle between 15° and 
164° and a transverse momentum exceeding 0.15 GeV, to reject non-e'^p background; 

• I Ay? I < 1 rad, where \Aip\ is the difference between the azimuth of the net transverse 
momentum as measured by the tracks associated with the vertex and the azimuth 
measured by the CAL. This requirement reduced the contamination from random 
coincidences of cosmic rays with e^p interactions; 

• PT,tracks/PT^^^ > 0-1) whcrc PT,tracks is the net transverse momentum of the tracks 
associated with the vertex (this condition was not applied if p™*^*^ > 25 GeV). This 
cut rejected events with additional energy deposits in the CAL not related to e^p 
interactions (mainly cosmic rays) and beam-related background in which p™"*** has a 
small polar angle; 

• events were removed from the sample if there was an isolated electron or positron 
candidate with energy above 10 GeV, to reject NC DIS events; 

• events were rejected if -Ebhac2/-Ebcal > 0.5 for -Ebcal > 2 GeV and -E^bhaci/-£^bcal > 
0.85 for -Ebcal > 8 GeV, where -E'bhaci(bhac2) is the energy deposited in the first 
(second) HAG section of BGAL and -Ebcal is the total energy deposition in BGAL. 
These requirements rejected beam-related background; 

• tracking requirements were not applied if the highest-E';^* jet (see below) in the event 
had 77-'°* > 2; in such a case, a tighter p™'^^ cut of 20 GeV was applied; 

• > 200 GeV^, to ensure high trigger efficiency; 

• y < 0.9, to avoid the degradation of the resolution in near ?/ ~ 1. 

The selected events were visually inspected and a few remaining non-e^p background 
events (1.2% of the final sample), mainly cosmic-ray and halo-muon events, were removed. 

Jets were identified in the pseudorapidity (77) - azimuth (0) plane of the laboratory frame 
using the kr cluster algorithm [25] in the longitudinally invariant inclusive mode [26]. 
This algorithm combines objects with a small relative distance dij, 

dij = min(Er,i, Erjf ■ {{rji - rjjf + (0^ - 0^)^), 

where i?T,j, Vi ^^^d 0* the transverse energy, pseudorapidity and azimuth of object i. 
The axis of the jet was defined according to the Snowmass convention [27], where 77-'^* 
((/^i*^*) is the transverse-energy weighted mean pseudorapidity (azimuth) of all the objects 
belonging to that jet. The jets were reconstructed using the CAL and corrected [24,28] 
for detector effects to yield jets of hadrons. Events with at least one jet in the range 
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— 1 < r/j'^* < 2.5 were retained. The inclusive-jet sample contained N'^'^^ = 5335 (870) and 
iVP°^ = 2122 (2284) jets with E^.^* > 14 GeV in the e'p (e+p) data, where iV"^g(p°«) is the 
number of jets selected in events with negatively-(positively-)polarized lepton beams. The 
dijet sample was selected requiring the jets with the highest and second-highest -E^* to 
have E^r^^^ > 14 and E^^^'^ > 5 GeV, respectively. The three-jet sample was selected from 
the dijet sample by requiring the third-highest E}^^ jet to have E^r^^^ > 5 GeV. The e~p 
(e^p) dijet and three-jet samples contained 1117 (464) and 109 (30) events, respectively. 
Eleven events contained a fourth jet with E}^^ > 5 GeV. 



4 Monte Carlo simulation 

Samples of Monte Carlo (MC) events were generated to determine the response of the 
detector and to evaluate the correction factors necessary to obtain the hadron-level jet 
cross sections. The hadron level is defined by those hadrons with lifetime r > 10 ps. 
The generated events were passed through the Geant 3.21-based [29] ZEUS detector- 
and trigger-simulation programs [12]. They were reconstructed and analyzed by the same 
program chain as the data. The CC DIS events were generated using the Lepto 6.5 pro- 
gram [30] interfaced to Heracles 4.6.1 [31] via Djangoh 1.3 [32]. The Heracles pro- 
gram includes first-order electroweak radiative corrections, vertex and propagator terms, 
and two-boson exchange. The CTEQ5D [33] proton parton distribution functions (PDFs) 
were used. The QCD radiation was modeled with the color-dipole model (CDM) [34] by 
using the Ariadne 4.08 program [35] including the boson-gluon-fusion process. To study 
the systematic effect of the modeling of parton showers in the correction of the data, sam- 
ples of events were generated using the Lepto model which is based on first-order QCD 
matrix elements and parton showers (MEPS). For the generation of the Lepto-MEPS 
sample^, the option for soft-color interactions was switched off. In both cases, fragmen- 
tation into hadrons was performed using the Lund string model [36] as implemented in 
Jetset 7.4 [37,38]. 

The photoproduction background was estimated using resolved and direct samples gen- 
erated using the program Herwig 5.9 [39]. After all the selection cuts described in the 
previous section, the contribution from photoproduction events to the inclusive-jet sample 
was estimated to be smaller than 0.5% overall and amounted to ~ 2% in the lowest i?^* 
bin. The NC DIS contamination was estimated to be smaller than the photoproduction 
background. No background subtraction was performed. 

^ The program Lepto allows the generation of only QPM events when using the MEPS option for e'p 
collisions; thus, only positron samples were generated and used as an estimate of the systematic effect 
also for e~p collisions. 
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The jet search was performed on the MC events using the energy measured in the CAL 
cells in the same way as for the data. The MC samples provided a good description of 
the measured distributions of the kinematic and jet variables. The data distributions as 
functions of 77-''^*, -Ey*, and x for the inclusive-jet samples are shown in Figs.[2](e~p) and 
El {g^p) separately for negatively- and positively-polarized electron and positron beams. 
Figure m shows the 77-^°*^, t]-""*^, E}^^^ and E}^^"^ distributions for the dijet sample. These 
figures also show the MC simulations normalized to the number of jets in the data. 

The same jet algorithm was also applied to the final-state hadrons to obtain the predictions 
at the hadron level. To correct the data to the hadron level, multiplicative correction 
factors, defined as ratios of the measured quantities for jets (events) of hadrons over 
the same quantity for jets (events) at detector level, were estimated by using the CDM 
samples and applied to the inclusive-jet (dijet and three-jet) data distributions. The 
samples of Lepto-MEPS were used as an estimation of the uncertainty on the modeling 
of the parton shower. 

Parton-level predictions were also obtained by applying the jet algorithm to the MC- 
generated partons. These predictions were used to correct the fixed-order QCD calcula- 
tions (see Section [5]) to the hadron level. The program Heracles was used to correct the 
predicted cross sections to the electroweak Born level evaluated using the fine structure 
constant a = 1/137.035999, the Fermi coupling constant Gf = 1.1664 ■ 10~^ GeV~^ and 
the mass of the Z boson Mz = 91.1876 GeV to determine the electroweak parameters. 

5 Fixed-order QCD calculations 

Fixed-order QCD calculations were obtained using the program Mepjet [40], which em- 
ploys the phase-space slicing method [41]. This is the only available program providing 
fixed-order QCD calculations for jet production in CC DIS. The jet algorithm described 
in Section [3] was also applied to the partons in the events generated by Mepjet to com- 
pute the predictions for the jet cross sections. The calculations were performed in the 
MS scheme. The calculations are 0{as) {0{a1)) for inclusive-jet (dijet and three-jet) 
production; this means that for inclusive-jet and dijet cross sections, the predictions are 
NLO whereas those for three-jet cross sections are only LO. The number of flavors was 
set to five and the renormalization (/i/j) and factorization (fip) scales were chosen to be 
fJ'R = fJ'F = Q- The calculations were performed using the ZEUS-S [42] parameterizations 
of the proton PDFs. Alternative calculations were performed using the CTEQ6 [43] and 
MRST2001 [44] sets of proton PDFs. The cross sections were evaluated using the same 
values for a, Gp and Mz as in the MC simulations (SectionH]). The mass of the W boson 
was fixed to 80.40 GeV [45]. The strong coupling constant was calculated at two loops 
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with = 226 MeV, corresponding to as{Mz) = 0.118. 

Since the measurements correspond to jets of hadrons whereas the QCD calculations 
correspond to jet of partons, the predictions were corrected to the hadron level using MC 
simulations. The multiphcative correction factor (Chad) is defined as the ratio of the cross 
sections for jets of hadrons to the same quantity for jets of partons, estimated using the 
MC programs described in Section |H The ratios obtained with the CDM model were 
taken as the default corrections, whereas those from the Lepto-MEPS model were used 
as an estimation of the effect of the parton shower. The value of Chad typically differs 
from unity by less than 5%, 10% and 30% for the inclusive-jet, dijet and three-jet cross 
sections, respectively. 

The following sources of uncertainty in the theoretical predictions were considered: 

• the uncertainty on the NLO QCD calculations due to terms beyond NLO, estimated 
by varying /i/j between Q/2 and 2Q, was typically below ±2% for the inclusive-jet 
cross sections and below ±5% for the dijet cross sections. For the LO calculations 
of the three-jet cross sections the uncertainty was ~ ±30%. For the three-jet cross 
sections, this uncertainty is dominant. Thus, no other theoretical uncertainty was 
taken into account for the three-jet cross sections; 

• the uncertainty on the NLO QCD calculations due to those on the proton PDFs was 
estimated by repeating the calculations using 22 additional sets from the ZEUS-S anal- 
ysis, which takes into account the statistical and correlated systematic experimental 
uncertainties of each data set used in the determination of the proton PDFs. The 
resulting uncertainty in the inclusive-jet e~p {e~^p) cross sections was below ±2 (4)%, 
except in the high-i?;^*, high-Q^ and high-x regions where it reached ±4 (10)%. The 
resulting uncertainty in the dijet e~p {e~^p) cross sections was below ±5 (5)%, except 
in the high-E';^*, high-Q^ and high-x regions where it reached ±7 (15)%; 

• the uncertainty on the NLO QCD calculations due to that on as{Mz) was estimated by 
repeating the calculations using two additional sets of proton PDFs, for which different 
values of as{Mz) were assumed in the fits. The difference between the calculations 
using these various sets was scaled by a factor such as to reflect the uncertainty on 
the current world average of [46]. The resulting uncertainty in the cross sections 
was below ±1%; 

• the uncertainty from the modeling of the QCD cascade was estimated as the differ- 
ence between the hadronization corrections obtained using the Ariadne and Lepto- 
MEPS models. The resulting uncertainty on the inclusive-jet and dijet cross sections 
was typically below 1%; 

• the uncertainty of the calculations due to the value of fip was estimated by repeating 
the calculations with fip = Q/2 and 2Q. The effect was negligible. 
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The total theoretical uncertainty was obtained by adding the individual uncertainties 
listed above in quadrature. 

6 Systematic uncertainties 

The following sources of systematic uncertainty were considered for the measurements of 
the jet cross sections [24]; values in percentage of the effects for the integrated inclusive-jet, 
dijet and three-jet cross sections are shown in parentheses: 

• the uncertainty in the absolute energy scale of the jets was estimated by studying the 
differences between data and MC simulation in single-jet NC DIS events [24, 47] by 
comparing the transverse energy imbalance between the scattered electron or positron 
and the jet. The uncertainty was found to be ±1% for < 75 GeV and ±2% for 
^j?* > 75 GeV (0.5%, 2%, 4%); 

• the uncertainty in the reconstruction of the kinematic variables due to the uncertainty 
in the absolute energy scale of the CAL was estimated by varying p™^*^*^ and yjB as 
measured with the CAL by ±1% for pf'' < 75 GeV and tl% for p^^"" > 75 GeV in 
the MC samples (0.2%, 0.3%, 0.4%); 

• the differences in the results obtained by using either CDM or Lepto-MEPS to correct 
the data for detector effects were taken as systematic uncertainties due to the modeling 
of the parton shower (0.7%, 7%, 6%); 

• the selection cut of p^^^^ > 11 GeV was changed to 10 GeV and 12 GeV in data and 
MC events (0.2%, 0.1%, below 0.1%); 

• the uncertainty in the simulation of the vertex position was estimated by changing the 
selection cut to -24 < Z^tx < 22 cm in data and MC events (2%, 2%, 2%); 

• the uncertainty in the simulation of the trigger was estimated to be negligible by using 
two different trigger selections in data and MC events (see Section [3]). 

The experimental uncertainties are dominated by the statistical uncertainty of the data, 
except for the inclusive-jet differential cross sections at high and high Q^, where 
the uncertainty coming from the modeling of the parton shower is large. The systematic 
uncertainties not associated with the absolute energy scale were added in quadrature to 
the statistical uncertainties and are shown in the figures as error bars. The uncertainties 
due to the absolute energy scale of the jets and the CAL were added linearly, due to the 
large bin-to-bin correlation, and are shown separately in the tables. In addition, there 
was an overall normalization uncertainty of 3.5% from the luminosity determination and 
3 — 5% uncertainty on the polarization measurement. 



8 



7 Results 



7.1 Polarized inclusive-jet cross sections 

Differential inclusive-jet cross sections were measured in the kinematic regime > 200 
GeV^ and y < 0.9. The cross sections were determined for jets with E';^* > 14 GeV and 
-1 < T]^"^ < 2.5. 

The inclusive-jet differential cross sections as functions of rj^^^, E^t^, and x for negatively- 
and positively-polarized e~p {e~^p) collisions are shown in Fig. [5] ([6]) and Tables [T] to 
m The predictions of the NLO calculations are compared to the data in the figures. 
The lower parts of the figures show the ratio of the cross sections for negatively- and 
positively-polarized lepton beams, which is in agreement with the measured polarization 
ratio, (1 - P"5)/(l - PP°^) = 1.79 ± 0.05 for e-p and (1 + PjD/ll + P = 2.10lJi:?^ 
for e~^p. The integrated polarized inclusive-jet cross sections, cxjets, are shown in Table [5l 
The measured cross sections are in good agreement with the predictions of the SM as given 
by the NLO QCD calculations, also shown in Table [5l in the kinematic range studied. 



7.2 Unpolarized inclusive-jet cross sections 

Figure [7] and Tables [6] to [13] show the unpolarized inclusive-jet differential cross sections 
as functions X in CC e^p DIS. The unpolarized cross section for an 

observable A was obtained via 



da ^ Kr\e^) iV^ 
dA^ ' Ce±-AA N^^' 



where 



data / ~ II pnog 1 ppos ) 

Ce± is the total integrated luminosity for e^p collisions and AA is the bin width. 

The measured da/drf^^ has a maximum at rf^^ ^ 1. The measured da/dE^^^ exhibits a 
fall-off of two (three) orders of magnitude in the e~p {e'^p) sample. Values of E"^* of more 
than 100 GeV are accessible with the present statistics. For 200 < < 2000 GeV^, the 
distributions display a weak dependence on Q"^. The cross sections as functions of i?;^* 
and show a less rapid fall-off than what is observed in NC DIS processes due to the 
presence of the massive W propagator. Furthermore, the measured cross sections for the 
e'^p sample decrease more rapidly as a function of E^^^ and than for the e~p sample 
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(see below). The values in x accessible by the data are within the range 0.013 < x < 0.63, 
as shown in Fig.[71i. The NLO QCD predictions using the ZEUS-S PDF sets are compared 
to the data in Fig. [71 Figure [8] shows the relative difference between the data and the 
predictions. The NLO QCD predictions give a reasonable description of the shape and 
normalization of the data. 

Figure[8]also displays the ratio of the e~p and e~^p differential cross sections. The measured 
ratio as a function of rf'^^ is constant and 2 as predicted by QCD, since the two reactions 
probe a different parton content of the proton. The ratio as a function of E";^* (Q^) 
increases as i?;^* (Q^) increases, in good agreement with the prediction. The increase at 
high values of E';^* and is expected due to the increasing contribution from the valence- 
quark densities in the proton at high x and the fact that both reactions are sensitive to 
different quark flavors. The behavior observed in the ratio of the measured cross sections 
as a function of is similar to the ratio of u and d parton densities. The same behavior 
is observed ClS db function of x. 

Figure [9] shows the contributions to the theoretical uncertainty from the terms beyond 
NLO, the parton-shower model and that coming from the uncertainty in the PDFs sep- 
arately for e~p and e~^p collisions. Also shown are calculations using other PDF sets. 
For inclusive-jet e^p CC cross sections, the uncertainty coming from that on the PDFs 
is dominant. At high i?;^*, and x, the uncertainty in the predicted cross sections for 
positron beams is larger than those for electron beams. This difference in the uncertainty 
due to the PDFs in the calculations for e~ and e"*" beams can be attributed to the different 
flavor content probed: in e~p {e~^p) at high x the W~ {W~^) will couple predominantly to 
the u (d) valence quark in the proton; at present, the uncertainty in the d parton density 
is larger than that for the u quark. Furthermore, the comparison with the calculations us- 
ing other PDF sets shows a wide spread in the predictions, especially for positron beams. 
Therefore, these measurements, in a phase-space region where the other theoretical un- 
certainties are well under control, have the potential to constrain the flavor content of 
the proton if used together with other data in global fits. A fast and accurate method 
to perform fits to extract the proton PDFs on data sets that included jet cross sections 
in NC DIS and photoproduction was recently developed by the ZEUS Collaboration [48]; 
the result was a sizable reduction of the uncertainty on the gluon density at medium and 
high X. Using the data presented here and extending such a method to jet cross sections 
in CC DIS may help to constrain the u and d valence quark distributions at high x. 

The integrated unpolarized inclusive-jet cross sections, Ujets; are shown in Table UM The 
measured cross sections are in good agreement with the predictions of NLO QCD, also 
shown in the table using different PDF sets. 
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7.3 Dijet cross sections 



Unpolarized dijet differential cross sections were measured in the kinematic regime > 
200 GeV^ and y < 0.9. The cross sections were determined for jets with E^r^^^ > 14 GeV, 
^j?*^ > 5 GeV and -1 < r]^"^ < 2.5. Figure M and Tables to [IH [15] and [IS] show 
the unpolarized dijet differential cross sections as functions of rf'^^, , and the dijet 
invariant mass, M^J, where r/j^* = (r/j<^" + r/j^*2)/2 and = (E^^" + E}^^'^)/2 in CC e^p 
DIS. The measured r/j*^* cross section has a maximum at rf'^^ ~ 1.25. The measured cross 
section as a function of Erp exhibits a fall-off of two orders of mag nitude for E^ > 20 GeV. 
For 200 < ^ 2000 GeV^, the distribution displays a weak dependence on Q^. Values 
of M-ij above 100 GeV are accessible with the present statistics. 

The NLO QCD predictions are compared to the data in Fig. [10] Figure [11] shows the 
relative difference to the predictions. The NLO predictions do not give an adequate 
description in shape and normalization of the measured differential cross sections over 
the entire phase space considered. In particular, for M^^, the data tend to be above 
the predictions for M^j > 70 GeV. It was reported [49] that calculations of jet cross 
sections in NC DIS computed using the Mepjet program differ by 5 — 8% from the 
results from other NLO programs. Comparisons of inclusive-jet calculations for NC DIS 
in the kinematic range of the measurements presented here performed using Mepjet and 
DiSENT [50] showed an agreement better than 1%. However, similar comparisons for 
dijet cross sections showed relative differences above ~ 5%. For CC DIS, it is not possible 
to quantify the degree of accuracy of the calculations of Mepjet since no alternative 
program exists. The NLO predictions give a reasonable description of the ratios of the 
cross sections for e~p and e~^p interactions (see Fig. [TTj) . New implementations of the 
theory are crucially needed to use the differential dijet cross sections presented here in 
global fits to extract the proton PDFs. 

The integrated unpolarized dijet cross sections are shown in Table [H] The measured 
cross sections are larger than the predictions of NLO QCD. 

7.4 Measurements of three-jet cross sections and observation of 
four- jet events 

Differential three-jet cross sections were measured in the kinematic regime > 200 GeV^ 
and y < 0.9. The cross sections were determined for jets with E}^^^ > 14 GeV, i?^*^ > 5 
GeV, ^j?*^ > 5 GeV and -1 < < 2.5. Three-jet cross sections in CC DIS were 
measured for the first time in e^p collisions. 

Figure [12] shows a three-jet candidate event in the ZEUS detector: a clear three-jet 
topology and large transverse momentum are observed. The three-jet selected sample 
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also contains 9 e~p and 2 e^p candidates with a fourth jet of transverse energy above 5 
GeV in the 77-'''* range considered. One of these candidates is displayed in Fig. [131 the 
fourth jet is clearly observed in the ZEUS detector. 

Figure [TH and Tables [6] to [HI [15] and [16] show the unpolarized three-jet differential cross 
sections as functions of ff'^^, E'^^^ , and the three-jet invariant mass, where ff^^ = 

above 100 GeV are accessible with the present statistics. 

The predictions of LO QCD are compared to the data in Fig. [H] The currently available 
QCD calculations are only lowest order and cannot predict the normalization of the data. 
Therefore, they were scaled by 1.92 and 1.42 for e~p and e~^p collisions, respectively, so as 
to reproduce the measured integrated three-jet cross section. The scaled LO calculations 
give a good description of the shape of the data. Figure [TH] shows the relative difference 
between the data and the scaled predictions. The lower part of Fig. [15] shows the ratio of 
the differential cross sections for the e~p and e~^p samples. 

The integrated unpolarized three-jet cross sections are shown in Table UM The predictions 
of LO QCD are also shown in the table. 

8 Summary 

Measurements of polarized and unpolarized integrated and differential multi-jet cross 
sections in CC e^p DIS were made using 0.36 fb~^ of data collected with the ZEUS 
detector at HERA II. The measurements were made in the kinematic region defined by 

> 200 GeV^ and y < 0.9. Jets were identified in the laboratory frame using the 
cluster algorithm in the longitudinally invariant inclusive mode. 

Polarized inclusive-jet cross sections were measured integrated over the phase-space region 
considered and differentially as functions of rj^^^, E^t^, and x for jets with i?^* > 14 
GeV and —1 < rf'^^ < 2.5. The measured cross sections are in good agreement with 
the SM predictions. The ratios of the differential cross sections for negative and positive 
longitudinally-polarized lepton beams are also well described by the predictions. 

Unpolarized differential inclusive-jet cross sections were measured as functions of rf^^, 
Et^i and x. The ratio of the differential cross sections for e~p and e^p collisions as a 
function of r]^^^ is ~ 2 in the t]-"^* range measured, as predicted by the theory. The ratio as 
a function of E-!^'^ (Q^) increases as E^^^ (Q^) increases, in agreement with the expected 
increased contribution from the valence-quark densities in the proton at high x and the 
fact that both reactions are sensitive to different quark flavors. Dijet differential cross 
sections were measured for jets with E^^"^^ > 14 GeV, E^l^^'^ > 5 GeV and —1 < rf'^^ < 2.5. 
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Next-to-leading-order QCD predictions computed using the program Mepjet were com- 
pared to the data. The NLO QCD predictions give a good description of the measured 
inclusive-jet cross sections. A detailed study of the theoretical uncertainties was per- 
formed: they are dominated by the contribution from the PDFs. Furthermore, the un- 
certainties due to the PDFs are larger for e^p than for e~p collisions. Therefore, these 
measurements, if used together with other data in global fits, have the potential to con- 
strain the flavor content of the proton at high x. 

The comparison of the predictions with the measured dijet differential cross sections shows 
a poor agreement in shape and normalization. Improved implementations of the theory 
are crucially needed to use these dijet measurements in a global fit to constrain the proton 
PDFs. 

Three-jet differential cross sections were measured for the first time in e^p collisions for 
jets with ^j?" > 14 GeV, > 5 GeV, E}^^^ > 5 GeV and -1 < < 2.5. The 

leading-order QCD predictions give a good description of the shape of the data. The 
three-jet sample also contains eleven candidates with a fourth jet of i?;^* > 5 GeV in the 
^jet range considered. 
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19.41 


±0.87 


±2.36 


+0 1 5 

_1_ W . X KJ 




P 


= +0.30 - 


inclusive jets 




— 1.0, —0.5 




2.81 


±0.37 


±0.21 


+0.11 
-0.10 


—0.5, 0.0 




7 99 


±0.55 


±0.28 


±0.16 


0.0, 0.5 




1 9 45 

X Zj . rtij 


±0.67 


±0.47 


+0.12 
-0.11 


0.5, 1.0 




1 fi 05 


±0.74 


±0.48 


±0.05 


1.0, 1.5 




1 5 QQ 


±0.74 


±0.50 


±0.05 


1.5, 2.0 




1 5 34 


±0.75 


±0.42 


+0.06 
-0.05 


2.0, 2.5 




1 1 3Q 

X X . OC/ 


±0.80 


±1.37 


±0.09 




P , 

e+ 




inclusive jets 








1.26 


±0.28 


±0.14 


+0 04 


-0.5, 0.0 




3 33 


±0.38 


±0.16 


±0.05 


5 




5 fi8 

?J .uo 


±0.47 


±0.14 


+0 04 


5 10 




fi fi4 


±0.49 


±0.40 


+0 04 


10 15 




fi 74 


±0.49 


±0.24 


+0 03 


1.5, 2.0 




fi 52 


±0.50 


±0.07 


±0.04 


2.0 2.5 




5.09 


±0.55 


±0.93 


+0 06 

—1— W . \J\J 






= +0.32 - 


inclusive jets 




-1.0, -0.5 




2.92 


±0.36 


±0.30 


±0.09 


-0.5,0.0 




7.37 


±0.49 


±0.29 


±0.11 


0.0,0.5 




10.34 


±0.54 


±0.41 


±0.08 


0.5,1.0 




14.75 


±0.63 


±0.65 


±0.08 


1.0,1.5 




12.38 


±0.58 


±0.65 


±0.06 


1.5,2.0 




11.64 


±0.58 


±0.50 


+0.07 
-0.06 


2.0,2.5 




10.22 


±0.68 


±1.26 


±0.12 



Table 1: Differential 'polarized inclusive-jet cross- sections dajdrf^^ for jets of 
hadrons in the laboratory frame selected with the longitudinally invariant clus- 
ter algorithm. The statistical, uncorrelated systematic and energy-scale (E'^) un- 
certainties are shown separately. 
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pj^* bin (CyeV) 


da/dE}^' (pb/GeV) 4tat 


"syst 


"ES 




P — — 27 - 


i n p1 n t;i \t(^ \ pf 






14, 21 


2 318 


+0 076 

—l—VJ . VJ ( \J 


±0 140 

— 1— VJ . X^VJ 


+0.045 
-0.051 


21, 29 


1 776 

X.I 1 \J 


+0 056 

_1_VJ . VJtJVJ 


±0.079 


+0.019 
-0.018 


29, 41 


1.419 


+0 038 


±0.070 


+0.009 
-0.010 


41, 55 


8S4 


+0 097 

—l—VJ .yjjLji 


±0.051 


+0.011 
-0.010 


55, 71 


450 


+0 01 8 

— 1— VJ .VJ ± o 


±0.022 


±0.011 


71, 87 
87, 120 


0.212 

0491 


±0.012 
+0 0035 

-I— VJ .VJVJ (JO 


±0.057 
±0.0225 


+0.003 
-0.002 
+0.0063 
-0.0050 




P = +0 30 - 








14, 21 


1.284 


+0 068 

—l—VJ . VJ\J<J 


±0.089 


+0.025 
-0.028 


21, 29 


1 081 

± .UO X 


+0 059 

-LVJ .VJtJzlj 


±0.049 


+0.012 
-0.011 


29 41 


0.795 


±0.035 


±0 041 

— 1— VJ . X 


+0.005 
-0.006 


41,55 


486 


±0.024 


±0.030 


±0.006 


55, 71 
71, 87 


263 
1 06 


+0 01 7 

—l—VJ . VJ X 1 

+0 01 

—l—VJ . VJ X VJ 


±0.016 
±0.029 


+0.007 
-0.006 

±0.001 


87 120 


0976 


+0 0035 

—l—VJ . VJVJOtJ 


±0 0149 

—l—VJ . W X 


+0.0041 
-0.0033 




P + — — 37 - 


IT"! fl n t;i A^P i pf 






14, 21 


761 


+0 053 

—l—VJ . VJtJ(J 


±0 058 

— 1— VJ . VJ tJ(J 


±0.012 


21, 29 
29, 41 


584 

VJ . 

338 

VJ . OtJO 


±0 039 

—l—VJ . VJOiJ 

±0 023 

—l—VJ . VJi^O 


±0.031 
±0.019 


+0.003 
-0.004 
+0.004 
-0.003 


41, 55 


1 60 

VJ . X VJVJ 


+0 01 4 

—l—VJ . VJ Xt: 


±0.011 


±0.004 


55, 71 


0609 


+0 0079 

^VJ .VJVJ ( C/ 


±0 0037 


+0.0025 
-0.0023 


71 87 


0.0157 


±0.0039 


±0 0047 

—l—VJ . IJVJ^ 1 


+0.0007 
-0.0006 


87,120 


00166 

VJ . VJW X VJVJ 


±0 00083 

— 1— VJ . \J\J\J(JfJ 


±0.00106 


+0.00032 
-0.00025 




p + = +0 32 - 


inclusive jets 






14, 21 


1.576 


±0.066 


±0.112 


+0.026 
-0.024 


21,29 


1.077 


±0.045 


±0.020 


+0.006 
-0.007 


29,41 


0.675 


±0.028 


±0.037 


+0.007 
-0.006 


41,55 


0.307 


±0.017 


±0.020 


±0.007 


55,71 
71,87 


0.1135 
0.0374 


±0.0094 
±0.0052 


±0.0143 
±0.0083 


+0.0047 
-0.0043 
+0.0016 
-0.0014 


87, 120 


0.0050 


±0.0012 


±0.0027 


+0.0010 
-0.0008 



Table 2: Differential polarized inclusive-jet cross-sections da/ dEi^ . Other details 
as in the caption to TableUi 
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bin (GeV2) 


da/dQ^ (pb/GeV2) 


<^stat 








Pe-=- 


-0.27 


- inclusive jet 


S 




200,500 
500, 1000 
1000,2000 


0.0246 
0.02226 
0.01578 




±0.0013 
±0.00079 
±0.00046 


±0.0015 
±0.00090 
±0.00053 


-1-0.0024 
-0.0023 
-1-0.00061 
-0.00059 

±0.00022 


2000,4000 


0.00797 




±0.00023 


±0.00007 


±0.00004 


4000, 10000 


0.002532 




±0.000072 


±0.000134 


±0.000038 


10000, 20000 


0.000429 




±0.000022 


±0.000029 


±0.000018 


20000,88000 


1.4- 10-5 




±1.3 - 10-6 


±5.7-10-6 


±1.4 - 10-6 






-0.30 


- inclusive jet 


s 




200,500 


0.0135 




±0.0011 


±0.0011 


±0.0013 


500, 1000 


0.01332 




±0.00074 


±0.00051 


+0.00036 
-0.00035 


1000,2000 


0.00926 




±0.00042 


±0.00030 


±0.00013 


2000,4000 


0.00437 




±0.00020 


±0.00004 


±0.00002 


4000, 10000 


0.001364 




±0.000064 


±0.000075 


±0.000020 


10000,20000 


30.6-10-5 




±2.2 - 10-5 


±2.0 - 10-5 


±1.3 - 10-5 


20000,88000 


8.1 - 10-6 




±1.2 - 10-6 


±3.3 - 10-6 


±0.8 - 10-6 




Pe+ = - 


-0.37 


- inclusive jet 


s 




200,500 
500, 1000 


0.00984 
0.00814 




±0.00096 
±0.00058 


±0.00049 
±0.00026 


+0.00091 
-0.00087 

±0.00016 


1000,2000 


0.00463 




±0.00030 


±0.00022 


±0.00002 


2000,4000 


0.00182 




±0.00013 


±0.00007 


±0.00003 


4000, 10000 


34.1 - 10-5 




±3.1 - 10-5 


±2.0 - 10-5 


±1.4 - 10-5 


10000,20000 


25.0 - 10-6 




±5.9 - 10-6 


±5.7- 10-6 


±2.3 - 10-6 




= +0.32 


- inclusive jets 




200,500 


0.0198 




±0.0012 


±0.0010 


+0.0018 
-0.0017 


500, 1000 


0.01586 




±0.00070 


±0.00023 


+0.00030 
-0.00031 


1000,2000 


0.00905 




±0.00036 


±0.00027 


+0.00004 
-0.00004 


2000,4000 
4000, 10000 


0.00345 
68.8-10-5 




±0.00016 
±3.9 - 10-5 


±0.00014 
±5.6 - 10-5 


+0.00005 
-0.00005 

±2.7-10-5 


10000, 20000 


61.3 - 10-6 




±8.0 ■ 10-6 


±10.0 - 10-6 


±5.6 - 10-6 


20000,88000 


3.5- 10-^ 




±1.8 ■ 10-^ 


±2.8 - 10-^ 


±0.6 - 10-^ 



Table 3: Differential polarized inclusive-jet cross-sections da / dQ'^ . Other details 
as in the caption to TableUi 
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X bin 


da j d 


X (pb/GeV) 5stat 










= -0.27 - 


inclusive jet 


S 




0.006,0.025 
0.025,0.063 




16.17 
46.9 


±0.70 
±1.3 


±0.83 
±1.6 


+1.06 
— 1.01 

±0.9 


0.063,0.16 




60.6 


±1.4 


±0.9 


±0.2 


0.16,0.40 




39.7 


±1.1 


±1.7 


±0.7 


0.40, 1.0 




3.94 


±0.33 


±0.53 


±0.29 






= +0.30 - 


inclusive jets 




0.006,0.025 




9.71 


±0.65 


±0.76 


+0.64 
—0.61 


0.025,0.063 




24.3 


±1.1 


±0.9 


±0.5 


0.063,0.16 




36.6 


±1.3 


±0.6 


±0.1 


0.16,0.40 




23.2 


±1.0 


±1.2 


±0.4 


0.40, 1.0 




2.21 


±0.30 


±0.38 


±0.16 




Pe+ = 


= -0.37 - 


inclusive jet 


s 




0.006,0.025 
0.025,0.063 




6.10 
13.59 


±0.54 
±0.81 


±0.28 
±0.53 


+0.35 

—0.34 

±0.15 


0.063,0.16 




14.90 


±0.82 


±0.62 


±0.16 


0.16,0.40 




5.22 


±0.47 


±0.30 


±0.13 


0.40,1.0 




0.33 


±0.11 


±0.12 


±0.02 




Pe+ = 


= +0.32 - 


inclusive jets 




0.006,0.025 




12.34 


±0.67 


±0.56 


+0.71 

-0.69 


0.025,0.063 




28.1 


±1.0 


±1.0 


±0.3 


0.063,0.16 




26.57 


±0.95 


±0.88 


±0.29 


0.16,0.40 




11.55 


±0.60 


±0.65 


±0.28 


0.40,1.0 




0.50 


±0.12 


±0.06 


±0.04 



Table 4: Differential polarized inclusive-jet cross-sections dojdx. Other details 
as in the caption to TableUi 
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lepton and polarization 


CTjcts (pb) 


4tat (pb) 


4yst (pb) 


^Es (pb) 


SM prediction (pb) 


Pe- = -0.27 ±0.01 
Pe- = +0.29 ± 0.01 


70.54 
40.53 


0.97 
0.88 


0.58 
0.45 


+0.43 
-0.40 
+0.24 
-0.23 


69.17 
38.67 


e+ ~ ^■'^ ' -0.02 

Pe+ = +0.32 ± 0.01 


17.55 
34.51 


0.60 
0.72 


0.57 
1.05 


±0.11 

+0.23 
-0.22 


16.86 
35.33 



Table 5: Integrated polarized inclusive- jet cross-sections a jets for jets of hadrons 
in the laboratory frame selected with the longitudinally invariant k^p cluster algo- 
rithm. The statistical, uncorrelated systematic and energy-scale fES ) uncertainties 
are shown separately. The uncertainty coming from the luminosity measurement 
is not shown. The predictions of the Standard Model as given by the Mepjet 
calculation are shown in the last column. 
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da/dr]i^^ (pb) 




^syst 


^ES 


Cqed 


Chad 


unpolarized - inclusive jets 


-1.0, -0.5 


3.58 


±0.28 


±0.30 


+0.15 
—0.13 


0.97 


0.93 


-0.5,0.0 


10.38 


±0.43 


±0.35 


+0.24 
—0.22 


0.97 


0.98 


0.0,0.5 


18.06 


±0.53 


±0.34 


±0.17 


0.98 


0.99 


0.5, 1.0 


22.17 


±0.57 


±0.75 


±0.07 


0.97 


1.00 


1.0,1.5 


22.38 


±0.58 


±0.90 


±0.07 


0.97 


1.00 


1.5,2.0 


20.33 


±0.57 


±0.31 


+0.08 

— U.U 1 


0.96 


1.01 


2.0,2.5 


15.59 


±0.61 


±1.85 


±0.12 


0.96 


1.01 


TtJ*^^ 1^111 
/ / Ulll 


da/dff''' (pb) 


^stat 






'-^QED 


"-^had 




unpolarized 


Hit g 
LllJcLo 








—1.0, —0.5 


0.103 


±0.073 


±0.003 


+0.033 
-0.032 


0.80 


0.60 


-0.5,0.0 


1.48 


±0.26 


±0.52 


±0.09 


0.96 


0.81 


0.0,0.5 


3.63 


±0.33 


±0.88 


±0.12 


0.97 


0.89 


0.5,1.0 


5.68 


±0.37 


±0.43 


+0.14 
-0.13 


0.98 


0.91 


1.0,1.5 


6.43 


±0.34 


±0.07 


+0.12 
-0.11 


0.97 


0.92 


1.5,2.0 


3.77 


±0.24 


±0.14 


±0.07 


0.95 


0.92 


2.0,2.5 


0.58 


±0.10 


±0.05 


±0.01 


0.93 


0.88 




unpolarized - 


three jets 








0.0,0.5 


0.49 


±0.18 


±0.49 


+0.03 
-0.05 


0.95 


0.75 


0.5,1.0 


1.05 


±0.21 


±0.15 


+0.05 
-0.06 


0.93 


0.78 


1.0,1.5 


1.06 


±0.17 


±0.13 


+0.04 
-0.03 


0.99 


0.80 


1.5,2.5 


0.246 


±0.046 


±0.053 


+0.009 
-0.007 


0.99 


0.80 



Table 6: Differential unpolarized inclusive-jet, dijet and three-jet cross-sections 
da/drj^^^ and da/dff^^ in e~p collisions for jets of hadrons in the laboratory frame 
selected with the longitudinally invariant kx cluster algorithm. The statistical, un- 
corrected systematic and jet- energy- scale (ES) uncertainties are shown separately. 
The multiplicative corrections for QED radiative effects, Cqed; one? the corrections 
for hadronization effects, Chad? he applied to the parton-level NLO QCD calcu- 
lations, are shown in the last two columns. 



22 



rf^^ bin 


da/dr]^''^ (pb) 


^stat 


^syst 


^ES 


Cqed 


Chad 


unpolarized - inclusive jets 


-1.0,-0.5 


2.12 




±0.24 


±0.21 


+U.U ( 
-0.06 


0.92 


0.94 


-0.5,0.0 


5.45 




±0.34 


±0.23 


±0.08 


0.94 


0.98 


0.0,0.5 


8.34 




±0.39 


±0.22 


±0.06 


0.95 


0.99 


0.5, 1.0 


10.90 




±0.43 


±0.54 


±0.06 


0.95 


1.00 


1.0, 1.5 


9.95 




±0.42 


±0.42 


±0.05 


0.95 


1.00 


1.5,2.0 


9.48 




±0.42 


±0.26 


+0.06 
-0.05 


0.94 


1.01 


2.0, 2.5 


7.88 




±0.48 


±1.15 


±0.09 


0.93 


1.02 


rp^^ bin 


da/drf''^ (pb) 


^stat 


^syst 


^ES 


Cqed 


Chad 






unpolarized 


- dijets 








-0.5,0.0 


0.92 




±0.23 


±0.33 


+0.07 
-0.05 


0.94 


0.83 


0.0,0.5 


1.75 




±0.28 


±0.42 


+0.06 
-0.05 


0.94 


0.90 


0.5,1.0 


3.01 




±0.32 


±0.24 


+0.07 
-0.06 


0.95 


0.94 


1.0,1.5 


3.52 




±0.28 


±0.06 


+0.07 
-0.06 


0.95 


0.94 


1.5,2.0 


1.95 




±0.20 


±0.11 


±0.04 


0.94 


0.94 


2.0,2.5 


0.287 




±0.073 


±0.043 


+0.008 
-0.005 


0.93 


0.91 






unpolarized - 


three jets 








0.0,0.5 


0.079 




±0.079 


±0.079 


+0.006 
-0.003 


0.92 


0.72 


0.5,1.0 


0.68 




±0.23 


±0.09 


+0.04 
-0.03 


0.93 


0.79 


1.0,1.5 


0.254 




±0.086 


±0.033 


+0.011 
-0.010 


0.94 


0.81 


1.5,2.5 


0.087 




±0.029 


±0.019 


±0.003 


0.91 


0.83 



Table 7: Differential unpolarized inclusive-jet, dijet and three-jet cross-sections 
da/drj^^^ and da/ff^^ in e^p collisions. Other details as in the caption to Table{^ 
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E}^' bin (GeV) 


da/dE}^' (pb/GeV) 


^stat 


^syst 




Cqed 


Chad 


unpolarized - inclusive jets 


14,21 


1.819 


±0.053 


±0.120 


+0.035 
-0.040 


1.03 


0.99 


21,29 


1.449 


±0.040 


±0.040 


+0.016 
-0.014 


1.01 


1.00 


29,41 


1.118 


±0.027 


±0.042 


+0.007 
-0.008 


0.97 


1.00 


41,55 


0.668 


±0.019 


±0.039 


±0.008 


0.95 


1.00 


55,71 


0.361 


±0.013 


±0.031 


±0.009 


0.92 


0.99 


71,87 


0.1596 


±0.0082 


±0.0360 


+0.0020 
-0.0015 


0.88 


0.99 


87,120 


0.0355 


±0.0026 


±0.0190 


+0.0053 
-0.0042 


0.83 


0.98 


bin (GeV) 


da/d'^T (pb/GeV) 


^stat 


^syst 


^ES 


Cqed 


Chad 


unpolarized - dijets 


9.5,14 


0.345 


±0.032 


±0.017 


±0.004 


1.03 


0.88 


14,21 


0.502 


±0.028 


±0.077 


+0.007 
-0.004 


0.99 


0.91 


21,29 


0.375 


±0.021 


±0.030 


±0.010 


0.94 


0.92 


29,41 


0.156 


±0.011 


±0.005 


±0.007 


0.92 


0.90 


41,55 


0.0390 


±0.0049 


±0.0107 


+0.0037 
-0.0026 


0.88 


0.90 


55,71 


0.0070 


±0.0023 


±0.0007 


+0.0014 
-0.0016 


0.90 


0.95 


71,87 


0.0021 


±0.0021 


±0.0021 


+0.0014 
-0.0000 


1.00 


0.94 


unpolarized - three jets 


8,9.5 


0.0096 


±0.0096 


±0.0025 


+0.0016 
-0.0017 


1.02 


0.72 


9.5,14 


0.060 


±0.015 


±0.015 


±0.001 


1.02 


0.77 


14,21 


0.099 


±0.015 


±0.003 


±0.004 


0.96 


0.79 


21,29 


0.0481 


±0.0087 


±0.0036 


+0.0030 
-0.0027 


0.92 


0.78 


29,41 


0.0078 


±0.0029 


±0.0008 


+0.0007 
-0.0006 


0.86 


0.85 



Table 8: Differential unpolarized inclusive-jet, dijet and three-jet cross-sections 
da/dE^rf^ andda/dlf^ in e~p collisions. Other details as in the caption to Table{^ 
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^j."* bin (GeV) 


da/dE}^' (pb/GeV) 




^syst 




Cqed 


Chad 


unpolarized - inclusive jets 


14,21 


1.200 




±0.046 


±0.087 


+0.020 
-0.019 


0.99 


1.00 


21,29 


0.863 




±0.033 


±0.029 


+0.005 
-0.006 


0.96 


1.00 


29,41 


0.522 




±0.019 


±0.025 


±0.005 


0.93 


1.00 


41,55 


0.242 




±0.012 


±0.016 


±0.006 


0.91 


0.99 


55,71 


0.0901 




±0.0067 


±0.0078 


+0.0037 
-0.0034 


0.87 


0.99 


71,87 


0.0269 




±0.0035 


±0.0067 


+0.0011 
-0.0010 


0.84 


0.98 


87, 120 


0.00329 




±0.00078 


±0.00179 


+0.00064 
-0.00050 


0.79 


0.97 


E^T bin (GeV) 


da/dEp (pb/GeV) 


<^stat 




^ES 


Cqed 


Chad 






unpol 


arized - dijets 








9.5,14 


0.288 




±0.033 


±0.017 


+0.001 
-0.002 


0.99 


0.90 


14,21 


0.330 




±0.026 


±0.052 


±0.004 


0.95 


0.93 


21,29 


0.190 




±0.018 


±0.015 


+0.007 
-0.006 


0.93 


0.93 


29,41 


0.0433 




±0.0066 


±0.0023 


+0.0026 
-0.0022 


0.88 


0.93 


41,55 


0.0112 




±0.0032 


±0.0031 


+0.0014 
-0.0010 


0.90 


0.95 


unpolarized - three jets 


8,9.5 


0.019 




±0.019 


±0.005 


±0.002 


0.82 


0.73 


9.5,14 


0.040 




±0.013 


±0.010 


±0.001 


0.96 


0.76 


14,21 


0.043 




±0.012 


±0.002 


±0.002 


0.93 


0.80 


21,29 


0.0034 




±0.0024 


±0.0003 


+0.0003 
-0.0002 


0.86 


0.82 


29,41 


0.0018 




±0.0018 


±0.0002 


+0.0003 
-0.0002 


0.88 


0.86 



Table 9: Differential unpolarized inclusive-jet, dijet and three-jet cross-sections 
da/dE^r^^ and da / dEf^ in e'^p collisions. Other details as in the caption to Tahle\^ 
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g2 bin (GeV^) 


da/dQ^ (pb/GeV2 


) ^stat 


^syst 


5es 


Cqed 


Chad 


unpolarized - inclusive jets 


200, 500 


0.01921 


±0.00087 


±0.00131 


+0.00190 
-0.00180 


0.98 


0.97 


500, 1000 
1000,2000 


0.01803 
0.01268 


±0.00056 
±0.00032 


±0.00047 
±0.00017 


+0.00049 
-0.00048 

±0.00018 


0.99 
0.98 


1.00 
1.00 


2000,4000 


0.00623 


±0.00016 


±0.00006 


±0.00003 


0.97 


1.00 


4000, 10000 


0.001963 


±0.000050 


±0.000097 


±0.000029 


0.95 


1.00 


10000,20000 
20000,88000 


0.000376 
1.121 ■ 10"^ 


±0.000016 
±0.091 ■ 10"^ 


±0.000045 
±0.706-10-5 


+0.000015 
-0.000016 

+0.113 in-5 
-0.114 ■ 


0.94 
0.93 


1.00 
1.00 






unpolarized - 


dijets 








200, 500 


0.00409 


±0.00091 


±0.00112 


+0.00094 
-0.00083 


1.00 


0.91 


500, 1000 
1000,2000 


0.00323 
0.00251 


±0.00031 
±0.00016 


±0.00037 
±0.00004 


+0.00027 
-0.00026 

±0.00009 


0.97 
0.97 


0.92 
0.91 


2000,4000 


0.001170 


±0.000071 


±0.000112 


+0.000024 
-0.000023 


0.96 


0.91 


4000, 10000 


0.000343 


±0.000021 


±0.000029 


±0.000011 


0.95 


0.89 


10000, 20000 


69.8 ■ 10"^ 


±6.9 ■ 10"^ 


±23.9 ■ lO-*' 


±4.1 ■ 10-^ 


0.99 


0.87 


20000,88000 


2.11 ■ 10-^ 


±0.37-10-^ 


±1.62 • 10-'5 


±0.26 ■ 10-5 


0.96 


0.86 


unpolarized - three jets 


500, 1000 


0.00034 


±0.00018 


±0.00015 


+0.00006 
-0.00005 


0.98 


0.78 


1000,2000 


0.000306 


±0.000073 


±0.000025 


+0.000027 
-0.000026 


0.96 


0.78 


2000,4000 


0.000193 


±0.000034 


±0.000024 


+0.000011 
-0.000010 


0.95 


0.79 


4000, 10000 


5.04-10-^ 


±0.93-10-5 


±2.18 ■ 10-5 


±0.20 ■ 10-5 


0.95 


0.79 


10000, 20000 


1.04-10-^ 


±0.25 ■ 10-5 


±0.70 ■ 10-5 


+0.09 1 n-5 
-0.08 ■ 


0.95 


0.75 


20000,88000 


1.34- 10-^ 


±0.99 ■ 10"^ 


±0.06-10-^ 


±0.19 • 10"^ 


0.92 


0.77 



Table 10: Differential unpolarized inclusive- jet, dijet and three-jet cross-sections 
da/dQ"^ in e^p collisions. Other details as in the caption to Tablel^ 
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bin (GeV2) 


da/dQ'^ (pb/GeV2 


) ^sta.t 






Cqed 


Chad 


unpolarized - inclusive jets 


200, 500 


0.01525 


±0.00082 


±0.00074 


-0.00134 


0.96 


0.97 


500, 1000 


0.01241 


±0.00050 


±0.00018 


±0.00024 


0.96 


1.00 


1000, 2000 


0.00707 


±0.00026 


±0.00021 


±0.00003 


0.95 


1.00 


2000, 4000 


0.00273 


±0.00011 


±0.00010 


±0.00004 


0.94 


1.00 


4000, 10000 


0.000530 


±0.000027 


±0.000037 


±0.000021 


0.91 


1.00 


10000, 20000 


4.35 • 10-5 


±0.53 • 10-5 


±0.71 • 10-5 


±0.40 - 10-5 


0.88 


1.00 


20000, 88000 


2.3 ■ lO"'^ 


±1.1 ■ lo-'^ 


±1.5 - 10-^ 


±0.4 - 10-^ 


0.81 


1.00 






unpolarized - 


dijets 








200, 500 


0.00204 


±0.00072 


±0.00051 


+0.00042 
-0.00037 


0.98 


0.92 


500, 1000 


0.00301 


±0.00033 


±0.00033 


+0.00025 
-0.00024 


0.96 


0.93 


1000, 2000 


0.00152 


±0.00014 


±0.00003 


+0.00004 
-0.00003 


0.95 


0.93 


2000, 4000 


0.000605 


±0.000058 


±0.000060 


±0.000016 


0.94 


0.92 


4000, 10000 


0.000106 


±0.000013 


±0.000011 


±0.000006 


0.91 


0.92 


10000, 20000 


9.3- 10-6 


±2.4-10-6 


±3.2 - 10-5 


+1.1 . in-6 
-1.0 


0.92 


0.92 


20000, 88000 


7.4 • 10"^ 


±7.4-10-^ 


±5.8 • 10-*^ 


±1.1 • 10-^ 


0.93 


0.88 


unpolarized - three jets 


500, 1000 


0.00038 


±0.00020 


±0.00017 


±0.00006 


0.97 


0.78 


1000, 2000 


0.000156 


±0.000055 


±0.000013 


+0.000014 
-0.000013 


0.91 


0.80 


2000, 4000 


0.000065 


±0.000021 


±0.000008 


+0.000003 
-0.000002 


0.91 


0.80 


4000, 10000 


9.2 ■ 10-6 


±4.4-10-6 


±4.1 - 10-6 


±0.7 - 10-6 


0.88 


0.79 


10000, 20000 


4.6 ■ lO-'^ 


±4.6-10-^ 


±3.1 ■ 10-^ 


±0.6 ■ 10-^ 


1.00 


0.82 



Table 11: Differential unpolarized inclusive- jet, dijet and three-jet cross-sections 
da/dQ^ in e^p collisions. Other details as in the caption to Tahle\^ 
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X bin 


da/dx (pb) 


^stat 


'^syst 


^ES 


Cqed 


Chad 


unpolarized - inclusive jets 


0.006,0.025 


13.12 


±0.50 


±0.74 


+0.86 
-0.82 


0.95 


1.01 


0.025,0.063 


35.82 


±0.85 


±1.13 


+0.70 
-0.69 


0.98 


1.00 


0.063,0.16 


49.26 


±0.98 


±0.09 


±0.14 


0.97 


0.99 


0.16,0.40 


31.81 


±0.77 


±1.30 


+0.53 
-0.54 


0.97 


0.99 


0.40,1.0 


3.11 


±0.23 


±0.30 


±0.23 


0.95 


0.98 



Table 12: Differential unpolarized inclusive-jet cross-sections da/dx in e p 
collisions. Other details as in the caption to Table{Bi 



X bin 


da/dx (pb) 


^stat 


^syst 




Cqed 


Chad 




unpolarized - 


inclusive jets 






0.006,0.025 


9.49 


±0.47 


±0.41 


+0.54 
-0.53 


0.93 


1.01 


0.025,0.063 


21.40 


±0.71 


±0.76 


±0.24 


0.96 


1.00 


0.063,0.16 


21.64 


±0.69 


±0.76 


+0.24 
-0.23 


0.95 


0.99 


0.16,0.40 


8.55 


±0.41 


±0.48 


±0.21 


0.93 


0.98 


0.40,1.0 


0.439 


±0.090 


±0.084 


±0.032 


0.91 


0.98 



Table 13: Differential unpolarized inclusive-jet cross-sections da/dx in e~^p 
collisions. Other details as in the caption to Table\B[ 
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lepton/ 


CTjets (pb) 


4tat (pb) 


5syst (pb) 


Ses (pb) 


QCD predictions (pb) 


jet multiplicity 










ZEUS-S 


CTEQ6 


MRST 


e~ /inclusive jet 


56.18 


0.68 


0.53 


+0.34 
-0.32 


54.47 ±0.75 


54.05 


54.56 


e"*" /inclusive jet 


26.88 


0.51 


0.82 


+0.18 
-0.17 


26.77 ±0.45 


25.85 


26.49 


e~ /dijet 


10.87 


0.34 


0.80 


+0.24 
-0.23 


9.14 ±0.35 


9.05 


9.26 


e"*" /dijet 


5.83 


0.29 


0.45 


+0.13 
-0.12 


4.57 ±0.19 


4.38 


4.55 


e~ /three jet 


1.52 


0.15 


0.09 


±0.06 


0.79 ±0.22 


0.79 


0.82 


e"*" /three jet 


0.563 


0.110 


0.037 


+0.025 
-0.022 


0.397 ±0.118 


0.386 


0.409 



Table 14: Integrated unpolarized jet cross-sections ajets for jets of hadrons in the 
laboratory frame selected with the longitudinally invariant kx cluster algorithm. The 
statistical, uncorrected systematic and energy-scale (ES ) uncertainties are shown 
separately. The predictions of QCD as given by the Mepjet calculations using the 
ZEUS-S PDFs are shown at NLO for the inclusive-jet and dijet cross sections and at 
LO for the three-jet cross sections, together with the total theoretical uncertainty. 
Also shown are the total cross sections predicted by QCD using the CTEQ6 or 
MRST PDF sets. 
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MJj bin (GeV) 


rfa/rfAfjj (pb/GeV) 


^stat 


^syst 




Cqed 


Chad 


unpolarized - dijets 


10,15 


0.090 


±0.011 


±0.010 


±0.001 


1.02 


0.87 


15,20 


0.207 


±0.018 


±0.011 


±0.002 


0.99 


0.88 


20,30 


0.284 


±0.016 


±0.032 


+0.005 
-0.004 


0.97 


0.89 


30,45 


0.250 


±0.014 


±0.028 


+0.006 
-0.007 


0.96 


0.90 


45,65 


0.0960 


±0.0083 


±0.0191 


+0.0051 
-0.0038 


0.95 


0.94 


65,90 


0.0385 


±0.0060 


±0.0058 


+0.0023 
-0.0032 


0.94 


0.96 


90,120 


0.0132 


±0.0046 


±0.0035 


+0.0023 
-0.0016 


0.97 


0.96 


M^j bin (GeV) 


da/dM'^^ (pb/GeV) 


^stat 


^syst 


^ES 


Cqed 


Chad 


unpolarized - three jets 


20,30 


0.0069 


±0.0024 


±0.0017 


±0.0001 


0.99 


0.83 


30,45 


0.0253 


±0.0045 


±0.0009 


+0.0006 
-0.0008 


0.99 


0.76 


45,65 


0.0307 


±0.0049 


±0.0012 


+0.0012 
-0.0011 


0.97 


0.77 


65,90 


0.0150 


±0.0034 


±0.0054 


+0.0014 
-0.0012 


0.94 


0.79 


90,120 


0.0036 


±0.0022 


±0.0001 


+0.0004 
-0.0003 


0.94 


0.85 



Table 15: Differential unpolarized dijet and three-jet cross-sections da/dM^^ and 
da/dM^^ in e~p collisions. Other details as in the caption to Tablel^ 
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Mjj bin (GeV) 


da/dMJj (pb/GeV) 


^stat 


^syst 


^ES 


Cqed 


Chad 


unpolarized - dijets 


10,15 


0.069 


±0.011 


±0.007 


±0.001 


0.96 


0.90 


15,20 


0.146 


±0.017 


±0.010 


+0.001 
-0.002 


0.96 


0.90 


20,30 


0.171 


±0.015 


±0.022 


±0.003 


0.95 


0.91 


30,45 


0.108 


±0.010 


±0.013 


±0.003 


0.95 


0.93 


45,65 


0.0522 


±0.0077 


±0.0103 


+0.0023 
-0.0025 


0.94 


0.95 


65,90 


0.0169 


±0.0052 


±0.0026 


+0.0015 
-0.0011 


0.94 


0.97 


90, 120 


0.0020 


±0.0020 


±0.0002 


+0.0003 
-0.0002 


0.92 


0.99 


M^j bin (GeV) 


da/dM^^ (pb/GeV) 


^stat 


^syst 


^ES 


Cqed 


Chad 


unpolarized - three jets 


20,30 


0.0051 


±0.0024 


±0.0013 


±0.0001 


0.91 


0.84 


30,45 


0.0124 


±0.0037 


±0.0006 


+0.0005 
-0.0003 


0.95 


0.75 


45,65 


0.0076 


±0.0029 


±0.0007 


±0.0004 


0.94 


0.78 


65,90 


0.0049 


±0.0024 


±0.0018 


+0.0005 
-0.0004 


0.90 


0.81 



Table 16: Differential unpolarized dijet and three-jet cross-sections da/dM^^ and 
da/dM^^ in e^p collisions. Other details as in the caption to Tahle\^ 
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(b) (c) 

Figure 1: Examples of jet production up to 0{as) in CC DIS. Feynman dia- 
grams for: (a) quark-parton model, (b) ho s on- gluon fusion and (c) QCD-Compton 
processes. 
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Figure 2: Detector-level data distributions for inclusive-jet production with neg- 
ative ( dots ) and positive ( open circles ) longitudinally polarized electron beams with 
E^r^^ > 14 GeV and —1 < rf^^ < 2.5 in the kinematic region given by > 200 
GeV^ and y < 0.9 as functions of (a) rf^^, (b) E}^^ , (c) and (d) Bjorken x. For 
comparison, the distributions of the CDM Monte Carlo model normalized to the 
number of jets in the data (solid histograms) are included. 
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Figure 3: Detector-level data distributions for inclusive-jet production with neg- 
ative ( dots ) and positive ( open circles ) longitudinally polarized positron beams with 
E^^ > U GeV and -1 < rj^^^ < 2.5 in the kinematic region given by > 200 
GeV^ and y < 0.9 as functions of (a) rj^^^, (b) E^r^^, (c) and (d) Bjorken x. The 
distributions of the Lepto-MEPS Monte Carlo model are also included (dashed 
histograms). Other details as in the caption to Fig. [3 
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Figure 4: Detector-level data distributions for dijet production with electron 
(dots) and positron (open circles) beams with E^^^^ > 14 GeV, E^r^^'^ > 5 GeV and 
— 1 < 7]^'^^ < 2.5 in the kinematic region given by > 200 GeV^ and y < 0.9 
as functions of (a) rj^'^^^ , (b) rj^^^'^ , (c) i?^*^ and (d) E^^^"^ . Other details as in the 
caption to Fig. [H 
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ZEUS 




Figure 5: Measured inclusive- jet cross sections in CC e p DIS for jets with 
E^r^^ > 14 GeV and —1 < 77-''^* < 2.5 in the kinematic regime given by > 
200 GeV^ and y < 0.9 as functions of (a) Tf""^, (h) E)^\ (c) and (d) Bjorken 
X for negative (dots) and positive (open circles) longitudinally polarized electron 
beams. The data points are plotted at the bin centres. The inner error bars represent 
the statistical uncertainties of the data, and the outer error bars show the statistical 
and uncorrected systematic uncertainties added in quadrature. For comparison, the 
predictions of NLO QCD (solid lines) are included. The lower parts of the figures 
display the ratio of the cross sections for negatively- and positively-polarized electron 
beams (squares) and the prediction (solid line) for the measured polarizations. The 
hatched band displays the uncertainty due to the polarization measurement. 
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Figure 6: Measured inclusive-jet cross sections in CC e^p DIS for jets with E^^^ > 



14 CeV and —1 < rf°^ < 2.5 in the kinematic regime given by > 200 GeV^ and 
y < 0.9 as functions of (a) rj^'^^, (b) 



and (d) Bjorken x for negative 



(dots) and positive (open circles) longitudinally polarized electron beams, 
details as in the caption to Fig. O 
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Figure 7: Measured unpolarized inclusive-jet cross sections in CC DIS for jets 
with E}^' > U GeV and -1 < rj^^^ < 2.5 in the kinematic regime given by > 
200 GeV^ and y < 0.9 as functions of (a) rj^^^, (h) E^r^^ , (c) Q"^ and (d) Bjorken 
X in e~p (dots) and e^p (open circles) collisions. For comparison, the predictions 
of NLO QCD based on the Mepjet calculations using the ZEUS PDF sets (solid 
lines) are also shown. Other details as in the caption to Fig. O 
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Figure 8: Relative difference between the measured cross sections of Fig.^and the 
corresponding NLO QCD predictions in e~p (dots) and e^p (open circles) collisions 
as functions of (a) ■rf'^^, (h) E^^^ , (c) Q"^ and (d) Bjorken x. The lower parts of the 
figures display the ratio of the cross sections for e~p and e'^p collisions (squares). 
The hatched areas display the theoretical uncertainty and the shaded areas display 
the uncertainty due to the absolute energy scale. Other details as in the caption to 
Fig. IB 
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Figure 9: Overview of theoretical uncertainties for the inclusive-jet cross sections 
in CC DIS for jets with E^^^ > lA GeV and —1 < 77^'^* < 2.5 in the kinematic regime 
given by Q'^ > 200 GeV^ and y < 0.9 as functions of (a) rf''^, (h) E}^\ (c) Q"^ and 
(d) Bjorken x in e~p (shaded areas) and e~^p (hatched areas) collisions. Shown are 
the relative uncertainties induced by the variation of the renormalization scale fiji, 
the uncertainties on the proton PDFs and hadronisation model. Also shown are the 
relative differences between the NLO QCD calculations using the CTEQ6 (dashed 
lines) or MRST (dot-dashed lines) PDF sets to the calculations based on the ZEUS 
sets. 
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Figure 10: Measured unpolarized dijet cross sections in CC DIS for jets with 
E^"" > 14 GeV, >h GeV and -1 < 77^'^* < 2.5 in the kinematic regime given 

by g2 > 200 GeV^ and y < 0.9 as functions of (a) rjj^S (b) E^^\ (c) and (d) 
M-ij in e~p (dots) and e'^p (open circles) collisions. Other details as in the caption 
to Fig. 
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Figure 11: Relative difference between the measured cross sections of Fig. [7 
and the corresponding NLO QCD predictions in e~p (dots) and e~^p (open circles) 

collisions as functions of (a) ff^^, (b) E^^ , (c) Q'^ and (d) M^j. The lower parts 
of the figures display the ratio of the cross sections for e~p and e^p collisions 
(squares). Other details as in the caption to Fig. O 
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Figure 12: Three-jet candidate event in CC DIS in the ZEUS detector. The 
energy deposition in the CAL is proportional to the size and density of shading in 
the CAL cells. The lego plots show the CAL transverse energy deposition projected 
in the rj — cf) plane. In the X — Y view, only the energy deposition in the barrel 
calorimeter is shown. 
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Figure 13: Four-jet candidate event in CC DIS in the ZEUS detector. The energy 
deposition in the CAL is proportional to the size and density of shading in the CAL 
cells. The lego plots show the CAL transverse energy deposition projected in the 
Tj — cj) plane. In the X — Y view, only the energy deposition in the barrel calorimeter 
is shown. 
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Figure 14: Measured unpolarized three-jet cross sections in CC DIS for jets with 
E^r^'^ > 14 GeV, E^r^'^ > 5 GeV, E^^^^ > 5 GeV and -1 < r/j^* < 2.5 in the 
kinematic regime given by > 200 CeV^ and y < 0.9 as functions of (a) fj^'^^ , 
(h) Et, (c) g2 and (d) M^J in e p (dots) and e^p (open circles) collisions. For 
comparison, the predictions of LO QCD based on the Mepjet calculations using 
the ZEUS PDF sets (solid lines) are also shown. The predicted cross sections are 
normalized to the measured three-jet cross sections (xl.92 for e~p and xl.42 for 
e~^p). Other details as in the caption to Fig. 
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Figure 15: Relative difference between the measured cross sections of Fig. 
and the corresponding LO QCD predictions in e~p (dots) and e+p (open circles) 

collisions as functions of (a) ff°^, (b) E^^ > (^) ^'^^ (^) ^'^^ ■ The lower parts 
of the figures display the ratio of the cross sections for e~p and e^p collisions 
(squares). Other details as in the caption to Fig. O 
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